1// noise from vortex-antivortex annihilation 
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The magnetic flux noise induced by vortices in thin superconducting films is studied. Rigid 
vortices in thin films as well as pancake vortices in the pancake gas regime are addressed. The 
vortex dynamics is described by a Feynman path integral which fully accounts for the balance 
between vortex entropy and vortex energetics. We find that vortex pair creation (annihilation) in 
the presence of pinning is a natural source for noise with a 1// spectrum. 
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Over the last decade it was noticed that the theory of 
strong interactions (QCD) and type II superconductors 
share a common feature: it was realized that the QCD 
vacuum sustains percolating vortices 0. The dynam- 
ics of the vortex texture is thereby related to important 
properties of QCD such as quark confinement. For in- 
stance, quark deconfinement at high temperatures ap- 
pears as vortex depercolation transition^, and the uni- 
versality class of the transition can be anticipated by in- 
vestigating the vortex dynamics only ^ . It turned out Q 
that a vortex model which essentially incorporates the 
vortex entropy and an energy penalty for vortex curva- 
ture already captures the dynamics of the confining vor- 
tices. Below, we will apply this idea to describe vortices 
of type II superconductors: In the case that the super- 
conductor partition function is dominated by the vortex 
entropy, the vortex dynamics is described by a Feynman 
"path integral" which correctly accounts for the vortex 
entropy. The only energy penalty is given for the vortex 
mobility. 

A central issue in the context of type II superconduc- 
tors is the magnetic flux noise and voltage noise gen- 
erated by the dynamics of vortices in superconductors. 
This topic has been investigated both theoretically and 
experimentally throughout the last four decades, start- 
ing with the pioneering work by Van Ooijen and Van 
Gurp IE • Early work focused on so-called flux- flow 
noise Q due to density and velocity fluctuations of vor- 
tices driven by a bias current, which causes voltage fluc- 
tuations in the flux flow regime with a noise power spec- 
trum often scaling with frequency / as 1//" (a « 1 — 2). 
More recently, various groups have performed experi- 
ments probing directly the magnetic flux noise in super- 
conductors which is sensed by nearby magnetic field or 
flux sensors 9, 10, IL,, 12, 13, Ijj. Within such an ap- 
proach, the magnetic flux noise at low magnetic fields 
has also been detected by placing a superconductor in 



close vicinity to a superconducting quantum interference 
device (SQUID) which senses directly the motion of vor- 
tices [13 • Similarly, one can simply detect the motion of 
vortices located in the SQUID structure itself 0|, e.g., 
combined with imaging of such vortices [T^ . Such mea- 
surements in low magnetic fields show also typically scal- 
ing of the spectral density of flux noise as 1//. 

Apart from the attempt to understand the origin and 
mechanisms of such noise as a prerequisite for improv- 
ing the properties of superconducting devices, the in- 
vestigation of noise was also driven by the wish to 
exploit vortex dynamics as a probe of the supercon- 
ducting state, in particular after the discovery of the 
high transition temperature Tc cuprate superconduc- 
tors, which show very rich phenomena associated both 
with the static and dynamic properties of 'vortex mat- 
ter' [TsL IT^ I20I l2l| . An interesting aspect is related 
to the layered structure of the cuprate superconductors 
and very thin superconducting films, which presumably 
leads to enhanced two-dimensional fluctuations near the 
Kosterlitz-Thouless transition as determined from sim- 
ulations of Xy-models with time-dependent Ginzburg- 
Landau dynamics [2^l23| . 

In this letter, the magnetic flux noise induced by a vortex 
to a slot, similar to the situation of vortices trapped in a 
SQUID washer [3, is studied on the basis of a statisti- 
cal vortex model. The model applies to the case of rigid 
vortices in thin films as well as to pancake vortices in 
the pancake gas regime close to Tc. Disregarding vortex 
pair creation (annihilation), we find that the noise spec- 
trum shows a l/^/J behavior. Taking into account vortex 
pair creation (annihilation), we observe a 1// law if the 
vortex pinning is active, while the spectrum behaves like 
1//^ if pinning is absent. This demonstrates that vortex 
pair creation (annihilation) of pinned vortices can be a 
natural source for 1// noise. 
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Let us consider a single "pancake" vortex which moves 
in the 2-dimensional xy plane. Vortex motion induces a 
time dependent magnetic flux to the slot parallel to the 
y-axis (see figure^ . The vortex "zitterbewegung" causes 
a frequency dependent noise to the SQUID, which mea- 
sures the flux induced to the slot. The derivation of the 
noise spectrum S{f) will be our goal below. The mo- 
tion of the vortex is described in the present case by a 
2-dimension vector x(t) which depends on the time t. In 
order to describe the pinning of the vortex to the ori- 
gin af = 0, we constrain the vortex motion by the weak 
pinning condition: 



i j'^ dt' x{t') 



0. 



(1) 



The latter condition implies that the "center of mass" 
of the vortex is tied to the origin of the coordinate sys- 
tem. The quantity of primary interest is the probability 
distribution 



P(xo,f;0,t) = (f(O)f(t)) 



(2) 



of finding a vortex at x at time t, when the vortex was 
located at xq at time t — 0. The average in ^ is taken 
over all vortex world lines which satisfy the pinning con- 
dition JQ) and the boundary conditions (xq , 0) and (x, t) . 
The noise spectrum is obtained from the probability dis- 
tribution by (w = 27r/) 



S{u>) = j dt (fxo (fx Td{x) P{xo, x;0,t) 



iut 



(3) 



Thereby, Td{x) is the transfer function which describes 
the coupling of the magnetic field of the vortex located 
a.t X = {x,y) at time t to the slot, which is loacted at dis- 
tance d from the y-axis. A simple model for the transfer 
function is provided by the choice 



Td[x) = 5(x - d) , 



(4) 



Thereby, it is assumed that only "pancake" vortices 
which possess overlap with the slot induce a signal. The 
X integration in |j2Jl sums up all contrubtions to the slot, 
whereas the xq integration can be understood as an av- 
erage over all initial conditions of the vortex. 

In the following, we will assume that the partition func- 
tion of the vortex matter is solely given by the inter- 
play between vortex entropy and vortex energy. In this 
regime, the dynamic of the vortices is captured by a sta- 
tistical model. Correlation functions are obtained by av- 
eraging over "pancake" vortex ensembles: The probabil- 
ity distribution is effectively described in terms of the 
"Feynman" path integral 

P(xo,f;0,t) - A/" y ^^(^^'^(7/ dt' x{t')^ (5) 



exp - / dt' -[i{t')] 



where m is the only parameter of the model. m 
parametrizes the mobility of the "pancake" vortex mov- 
ing in the 2d plane. All vortex worldlines, which con- 
tribute to the path integral |(SJ) start at xq at time t — 
and end at x at time t. The normalization constant M 
can be obtained from the condition J d^x P{xo, x; 0, t) = 
1 . Let us calculate the auxilliary quantity (the depen- 
dence on. Xq, X is suppressed) 



Z{a,t) = J Vx{t) expjzy^ dt' x{t'] 



(6) 



exp 



dt' 



Z(0,t) corresponds to the probability distribution of 
freely moving vortices, while integration over a installs 
the weak pinning condition (Q, i.e., 



P(xo,x;0, i) = 



da 



z{d,t) 



(7) 



It is tedious, but straightforward to calculate the function 
Z{d, t) in ((HJ. For this purpose, the time interval [0, t] is 
divided into N pieces of equal distance: 



t„ = — n 
N 



X{t'n) 



n = Q 

XN 



JV 



X 



The discretized version of the path integral jnj is given 
by 



Z(a,i) 



-N/2 



dxi 



N 

. . dxN-1 exp<{ i— ^ Xr, 

n=0 



N-l 



2t 



(^n-l-l - XnY 



n=0 



where the prefactor of the functional integral arises from 
the measure Vx(t). After performing the Gaussian inte- 
grals and taking the limit N 00, we obtain up to an 
unimportant numerical factor 



r 771 2"! 

Z[a,t) cx exp |-— (a; - xo) | 



(8) 



exp 



24777 



d^ + —d{x + Xq) 



Setting d = in the last expression, we find that Z only 
depends on x — xo refiecting translation invariance in the 
absence of the pinning center at the origin. Extending 
the correlation function to negative times, we replace t 
in |(SJ) by Inserting (jS)) in Q and introducing the 
relative momentum q which is the conjugate variable to 
x — Xq, the integration over t, xq and a can be performed 
yielding 



S{U!) CX 



d'^q 
(2^ 



d'^X^' 



1 



r_ 

6m 



Td(x) e 



(9) 
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For the idealized transfer function in we finally find 
for the case y/Gmuj d <C 1 that 

S{lu) oc \ — . (10) 

V UJ 

Hence, statistical ensembles, consisting out of single and 
weakly pinned vortices, cannot account for the 1// type 
spectrum. 

In order to allow for the merge of a vortex antivortex 
or for the vortex pair creation in the vortex statistical 
ensemble, the vortex world line is addressed by an implict 
parameterization: 

x{t) ^ (x{T),y{T)) , t{T) , x-~{x,y,t), 

where the parameter r will be called "proper time" be- 
low. If t{T) is a non-monotonic function, vortex pair 
creation and annihilation is take into account. The weak 
pinning condition ^ is generalized to an average over 
the vortex world line 



cLt' x{t') = 0. 



(11) 



The vortex statistical ensemble is now approached by 
means of the path integral 



ZcA{a,t) = j Vx(t) expji^^ dr' x{t')^ (12 



exp 



dr' 



—x^ + ^t^ 



where the dot denotes the derivative with t. If /i is small, 
t(r) possesses many non-monotonic parts. This implies 
that 1//1 is related to the vortex pair creation rate. The 
ensemble average must be taken over all vortex world 



lines starting at Xq :— (xo,?/o,0) and ending at x_ :— 
{x,y,t). There is no restriction for the proper time r 
implying that we will take the limit t oo below. A 
calculation, analogous to that in the previous subsection, 

yields 

f°° ( rn ^ ^ 2 M 2I 

ZcA{oi,t) cx j dr exp |-— (x - a;o) ~ | 



exp 



^—a^ + -d{x + x,)\. (13) 



After switching to the momentum space, the r integra- 
tion can be performed, i.e., 



ZcA{a,t) cx 



d'^q gni^-^o) giqotgi/2a(x+xo) 



(14) 



2m ' 2/i ' 24m 



The noise spectrum S{ljj) is now easily obtained by in- 
serting (fT^ in 0. Setting x = {x,y) and q — (91,(72), 
the a^o, <5, Qo and t integration is straightforward: 



S{uj) (X 



d'^q 



d^x Td{x) 



2q^ . 
3m 2fi 



(15) 



For the idealized transfer function in our final result 
is 



b{uj) (X exp<— W ujd> — 

V M J 2/^ 



(16) 



This is the central result of our letter: allowing for vor- 
tex pair creation and annihilation, the noise spectrum 
shows 1 /lu dependence at small frequencies if the vortex 
is pinned sufficiently close {utd <C 1) to the slot. 

How does pinning affect the noise spectrum? Let us con- 
sider the vortex statistical ensembles which include vor- 
tex pair creation and annihilation, but where the pinning 
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constraint is absent. The noise spectrum is given in this 
case by 

S{uj) oc J dt (fxo (fx Tdix) ZcAiO,t) e^"* . (17) 

Note that we have set a = in the auxihary function 
ZcA{ct,t) rather than integrating over a. Using (|14|l . 



the frequency dependent part of S{uj) is given by S{u}) oc 
l/w^. Obviously, the pinning of the "pancake" vortices 
is important to obtain a 1// noise spectrum. 
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